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Abstract 



o 
o 
o , 

CN \ The charge of the leading particle in DIS interactions is studied with the QCD models 

LEPTO, ARIADNE, RAPGAP and HERWIG. It is predicted to be preferentially positive, 
as is expected from the modified quark parton model and is studied over a wide range of 
t^. . x and Q 2 . 

m '■ 1 Introduction 

o . 

In deep-inelastic lepton-proton scattering processes, a small excess of positively charged lead- 
ing particles is expected owing to the valence quark structure of the proton. An interaction 
involving a valence u quark is 8 times more likely than that involving a valence d quark. No 
^ excess is expected from sea quark or gluon initiated interactions. 

Since leading hadrons "remember" the charge of the struck quark, this gives the experimentalist 
Ph! a probe to investigate the partonic structure of the proton. At HERA this method could be 
used to perform such studies over several orders of magnitude in x and Q 2 and will complement 
existing studies based on jet measurements and other hadronic final state variables [1] [2]. 



X. 

1.1 Definitions 

To define an event charge, let us first determine a reference direction on which to project mo- 
menta. If the DIS scattering were simply an elastic deflection of a quark by a positron, this 
direction would simply be the final quark direction. But we know that the reaction is more 
complicated, involving gluons or even more than one quark (boson-gluon-fusion BGF). Never- 
theless, from the deflected positron parameters (which, for simulated events are known without 
error) we define the 'hadronic axis' by assuming the scattering to be elastic, and using the 
deflected positron information. 

Alternative methods could be to use event shape information, such as the thrust axis, to define 
the hadronic system. This is beyond the scope of this paper and is not considered here. 
Event charge studies have already been made in deep inelastic scattering [3] and in e+e- inter- 
actions at LEP[4]. 

For this paper, the following definition of event charge is used: Q = J2i PiQi/Hi Pi ■, where 
is the particle charge (-1 or +1) and Pf is the projected momentum on the hadronic axis .The 
sums run over i all particles which have positive Pj"'s; Q is between -1 and +1. This was used 
by Aleph in [4]. 
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1.2 Relation of event charge Q to scattering on valence quarks. 
The Leading particle and the positive excess e + 

For this study, the leading particle is defined as the charged hadron which has the largest 
momentum projected onto the hadronic axis. If the experimentally measurable K® or A were 
to be the leading particle then the event is discarded. 

The positive excess e + is defined as the ratio between the number of leading particles which 
are positively signed iV + over those which are negatively signed N_: e + = j±. 

1.2.1 Relation of e + to interactions of valence quarks 

In this section, we will use e + to explain the relation of the positive mean event charge to 

interactions involving the proton valence quarks. The magnitude of e + can be related to two 
quantities: 

> the positive excess at the quark level, given by the ratio = -|p where M~ = + 
N d + N s and M+ = N U + N 1 + N-g; N u , N d , N s , N„, N 2 and Nj are the numbers of those 
events where the interaction has given a leading u-,d-,s-,TZ-,d-, respectively s-quark. 

D> p: the probability that the charge of the initial quark and the one of the leading particle are 
the same. Here p should be the same irrespective of the quark flavour. This will later be 
shown a flawed assumption, but we can take mean values. 

Then the dependence of e + upon e+ and p is described by the formula: 

e + = (pM + (l - p)M-)/{pM- + (1 - p)M + ) =(l-p + £+p)/(p + e+(l - p)). 

The probability p, averaged on different flavours, is given by QCD based models, and is 
found to be of the order of 0.6 to 0.7. 

The valence quark interactions from the proton are 8 times more likely to be initiated by an u 
quark than by a d quark . Hence, for valence quarks, more n + and K + than 7r~ and K~ should 
be detected as leading particles. 

For example, in the very high x region, the valence quark interactions will dominate the ep 
cross-section. Therefore, £+ is equal to 8, and e + is expected to be 1.39 for p = 0.6. 
To quantify this, we take Y+ and Y~ as the yields of positive and negative quarks struck by 
the photon, respectively, and produced by some process a. To simplify, we could consider only 
three different processes a: Boson Gluon fusion (BGF), scattering on a sea quark (Sea) or 
scattering on valence quarks (V). 

Here, we do not distinguish between elastic e — q diffusion and QCD Compton reactions, where 
one gluon is added: e + + q — > e /+ + q' + gluon. It has been shown with LEPTO, that so far 
as the event charge is concerned, they are equivalent. This can be simply understood: in QCD 
compton events, the gluon is often weaker than the final quark q', and also it has to break down 
into two quarks (plus gluons). The chance that one out of these two quarks has a larger energy 
than q' is very limited; so the leading particle will generally stay related to the q' quark. 
The positive excess at the quark level is thus: e+ = (Ybgf + Y s ta + ^v)/(Xbgf + Xsea + ^v)- 
For sea quarks and quarks emitted through BGF, the number of quarks is identical to the 
number of their antiquarks, we can state: 

YbgF + Ysta = YbGF + Y.sea = ^0- 
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Then: e+ = (1 + ^)/l + ^-). Let i? = ^ and 5 = \, so that ^- = |. As was quoted 

above, the value of R should be 8 (naively), so that : = (1 + S)/(l + |[) is nearly equal to 

1 + S. We see then that the positive excess at the quark level is a measure of the ratio of yield 
by valence quarks over yield by other processes. 

We access to it by measuring the mean event charge, but the probability p introduced above 
should be known, using QCD based models. We can perform some tests on these QCD based 
models, mainly on various fragmentation properties [5] and leading particle studies (mainly Kjjj 
and A [6]). After determining the mean event charge as a function of x and Q 2 , we cannot 
have direct values of e+ but we may compare the reconstructed charges with those predicted 
by QCD based models. 

1.3 Event selection 

The leading particle plays a major role in determination of the mean charge, and so we will 
perform our kinematic cuts on the leading particle parameters. In the following, we will cut 
all events having a leading particle angle less than 0.7 rad (where the angle is defined with 
respect to the initial proton direction ) or larger than 7r — 0.7 rad. In the forward direction, we 
then avoid particles coming from the proton remnant. In both directions, we have to reserve 
some room for the next to leading (n.t.l.) and other particles, so that the event charge is well 
evaluated. 

Another cut on all tracks has been imposed by experimental conditions at HI: to avoid beam 
gas background, to which the mean value of the event Charge Q is very sensitive, we require 
all tracks to have a transverse momentum larger than 0.6 GeV/c. 

The parameter / = P\\ ading / P qU ark "was also used to select events where the leading particle 
"remembers" better the initial quark flavour. Typically, it was required to be larger than 0.15. 

1.4 Our sample of events 

There exists different QCD based models, for which we take the following versions: LEPTO(6.5) [7], 
ARIADNE(410)[8], RAPGAP(2.0)[9], and HERWIG(5.9008)[10]. LEPTO, ARIADNE, and 
RAPGAP have the same hadronisation program, JETSET [11]. 

For each of these four QCD based models, we have generated 2 million events in the range 
3. < Q 2 < 100 GeV 2 . The corresponding luminosity is 6.3 pb _1 , which is below what is 
now available at HERA. The GRV [12] parton density was used for all. For ARIADNE and 
RAPGAP, the pomeron was off. For LEPTO, the QCD effects were on. 

2 Fragmentation studies 

We have made comparisons of different QCD based models for the following fragmentation 
parameters: 

- charge correlations in the system of the 2 or 3 leading particles; 

- projected momentum of the leading particle on the hadronic direction; 

- <f, rapidity and momentum ratio correlations of the leading and the struck quark; 

- ip, rapidity and charge correlations of the leading and n.t.l. particles; 

- Kg, A as leading particles; 

- p, K*^ and (f> production in the system of the leading and next to leading particles. 
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We show here charge correlations in the system of the 2 or 3 leading particles. However, 
the other parameters are shown in the Internet version of this note. 

All these quantities can be experimentally measured, and comparison to M.C's provides a good 
test of hadronisation programs. The detailed results are shown in the extended version (on the 
web). Generally, there is a good agreement between different QCD based models. 

2.1 Charge correlations for leading and n.t.l. particles 

The leading and next to leading particles may either have opposite charges, which is foreseen 
by the naive model of quark hadronisation, or same charges, which is also foreseen but in less 
frequent cases. In the following, when we write a pair of signs, the first is for the leading, the 
second for the n.t.l.. In the next figure 4 different ratios are shown, for channels 1 to 4: 



1: 



-and- 



++and- 



-■. this is the 'negative excess' of the n.t.l. Naively, this should be a little larger 
than 1, and it is very close to 1 for all QCD based models. 

2: this is the positive excess for leading and n.t.l. when both have the same charge. 

3: 3^: positive excess for opposite charge leading and n.t.l, lower than the preceding one. 

4: 333: this is the positive excess when the three first particles carry the same sign. 
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Figure 1: 4 ratios (defined in the text) for various x and Q 2 and for / = ' ead '" 9 > 0.15. 

r quark 
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The agreement between the three JETSET QCD based models and on the other hand HER- 
WIG is rather poor. 



In the next figure, we show 3 other ratios: 

1: dlfferent . slgns for the leading and the n.t.L.It should be much larger than 1, and it is of the 

same signs ° ° ' 

order of 2. 



2: 



3: 



different signs ^ leading and the next to next to leading (n.t.n.t.L). Naively, this would 

same signs ° ° v / J ' 

be expected to be less than 1 but is larger than 1. This is due to the fact that events 
with signs (+++) or ( ) are very seldom, as is shown in the next channel. 

(++-)+(+-+)+(+—)+(—+)+(-+-)+(-++) 

3*((+++)+( )) 

We see that the value of this ratio is between 3 to 4: it is very seldom that the three first 
particles carry the same sign. 
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Figure 2: 3 ratios (defined in the text just above)for some x and Q 2 and for / = > 0.15. 
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3 QCD model predictions for different parton induced 
quantities 



As was written in the introduction, the event charge Q depends upon two quantities: the mean 
probability p that the initial quark charge and that of the leading particle are the same and 
the positive excess at the quark level. Before presenting the results on the event charge, 
we will present the predictions of different QCD based models for these quantities. 



3.1 The sign excess depending on the flavour 



The sign excess related to a quark is defined as the ratio of the number of events originated 
from that quark, which gives a leading particle having its sign, to the number of those which 
originate from q as well, but give a leading particle of opposite sign. For instance, for the u 



quark initiating the reaction, it is: S u = whereas for the u:Su = jt- The sign excess 



N- 



N. 



N+- 



is related directly to the probability p: for example, if we have a sign excess of 3, then the 
probability p for that quark is 75 percent. Here we compare the S q 's for the 4 QCD based 
models. 



5.3 
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Sign excess as a function of the flavor of the struck quark 
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Figure 3: Sign excess as a function of the flavour of the struck quark, for some x and Q 2 . 

Each box of the preceding figure has 6 channels, from -3 to 3. We use the standard quark 
identity of the Particle Data Group [13]:-3 is s, -2 is u, -1 is d, 1 is d, 2 is u and 3 is s. The 
sign excesses are shown for these 6 abscissa of the boxes. 
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It is easily seen that the d,d have less sign excess than u,u; there are at least 2 reasons: first, 
in the naive model, the d can give a proton as leading particle,which then has the wrong sign; 
second, for what concerns charged particles, the u can give a n + or a K + , the d only a it~ and 
not a [14], and the u has there also an advantage. 

It is seen that the sign excesses for a quark and its antiparticle are equal (within statistical 
errors), which is a good consistency check of the models. 

Generally, HERWIG gives a greater sign excess than RAPGAP and LEPTO, which in turn 
gives more than ARIADNE. Let us point out that the sign excess reflects on one hand the 
hadronisation, (production of mesons and baryons), and on the other hand the production of 
gluons and quarks from an initial quark (the one from the proton) . As ARIADNE uses the 
same hadronisation program as RAPGAP and LEPTO, we see that there is a sensitivity to the 
underlying differences in the physics of these models. 

3.2 The positive excess at the level of initial quarks 

The positive excess at the level of quarks: e+ = (N u + N-^ + Nj) / (N^ + N d + N s ) , where N q is 
the number of events where a quark q from the proton initiates the reaction. On the next plot 
is shown how this quantity varies with x and Q 2 . We can see that, as expected, e q increases 
with x (more valence at higher x), and that, at large Q 2 , it begins to take important values. 
The different QCD based models show similar behaviours, except for HERWIG for 
Q 2 < 10 GeV 2 . 
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Figure 4: The positive excess at the level of initial quarks as a function of logi x and Q 2 
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3.3 Final result: mean event charge as a function of x and Q 2 



To recap, the charge that we consider here is the following: Q = 



y.p u 



, where the sum over % 



extends to all particles having a positive momentum projection on the struck quark direction. 
The found mean event charges are shown in the next figure. HERWIG predicts a far greater 
charge than the other models. The charge increases with x for the two ranges in Q 2 considered. 
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Figure 5: Event charge as a function of logiox and Q 2 for events within the cuts (mainly 
0.7 < leading particle 6 < n — 0.7). 



4 Conclusions 

Using four QCD based models the charge properties of inclusively produced hadrons in ep 
DIS processes has been investigated. Deviations of these models have been seen, particularly 
between those models employing string and cluster hadronisation for the total event charge. Es- 
pecially, HERWIG seems to give very different results than LEPTO, ARIADNE and RAPGAP, 
which all three use the hadronisation program JETSET. 
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The mean event charge may also be measured in special cases of deep inelastic scattering: 

-i- for charge current events, where it should have the sign of the incoming electron or 
positron in the beam. Due to limited statistics, this could only be an interesting check. 

-ii- for diffractive events. In the Pomeron model, the mean event charge is obviously zero. If 
we find a mean positive charge in these events, this means either that there is another 
mechanism involved in the diffraction or that the cuts selecting these diffractive events 
are not severe enough. 

Michel Jaffre and Ursula Berthon are acknowledged for their help in processing data. 
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